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Nonaxisymmetric Magnetorotational Instability in Proto-Neutron Stars 

Youhei Masada^’^, Takayoshi Sano^, and Hideaki Takabe^ 

ABSTRACT 

We investigate the stability of differentially rotating proto-neutron stars (PNSs) with a 
toroidal magnetic field. Stability criteria for nonaxisymmetric MHD instabilities are derived us¬ 
ing a local linear analysis. PNSs are expected to have much stronger radial shear in the rotation 
velocity compared to normal stars. We find that nonaxisymmetric magnetorotational instability 
(NMRI) with a large azimuthal wavenumber m is dominant over the kink mode {m = 1) in dif¬ 
ferentially rotating PNSs. The growth rate of the NMRI is of the order of the angular velocity fl 
which is faster than that of the kink-type instability by several orders of magnitude. The stability 
criteria are analogous to those of the axisymmetric magnetorotational instability with a poloidal 
field, although the effects of leptonic gradients are considered in our analysis. The NMRI can 
grow even in convectively stable layers if the wavevectors of unstable modes are parallel to the 
restoring force by the Brunt-Vaisala oscillation. The nonlinear evolution of NMRI could amplify 
the magnetic fields and drive MHD turbulence in PNSs, which may lead to enhancement of the 
neutrino luminosity. 

Subject headings: instabilities — MHD — stars : magnetic fields — stars : neutron 


1. INTRODUCTION 

Neutron stars are formed in the aftermath of supernova explosions associated to the gravitational core¬ 
collapse of massive stars (8 - 30Mq) at the end of their evolutions. The binding energies released by 
core-collapse are stored in the interior of a newly born hot neutron star, and those energies are emitted as 
the neutrinos in the explosion stage. Neutrino emission from proto-neutron stars (PNSs) is one of the key 
processes of the delayed explosion scenario of core-collapse supernovae, in other words, the size of neutrino 
luminosities and the energy deposition efficiency determine whether the explosion will succeed or not (Bethe 
1990). 

Convection in PNSs is considered to play an important role in enhancing the neutrino luminosities 
(Epstein 1979; Lattimer & Mazurek 1981; Arnett 1987; Burrows & Lattimer 1986; Burrow & Fryxell 1993). 
Other hydrodynamic instabilities may also contribute to the enhancement of neutrino luminosities (Bruenn & 
Dineva 1996; Mezzaacappa et al. 1998a; Miralles et al. 2000, 2002, 2004). Mayle (1985) and Wilson & Mayle 
(1988, 1993) have argued that the material in the PNSs would be unstable to a double diffusive instability 
which is referred to as “neutron finger”. Bruenn et al. (2004) have found a new double diffusive instability 
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called “lepto-entropy finger instability”. However, despite considering various hydrodynamic instabilities, 
recent numerical simulations have not succeeded in producing sufficient neutrino luminosities leading to the 
delayed supernova explosions (Burrows et al. 1993; Keil et al. 1996; Mezzacappa et al. 1998b; Buras et al. 
2003). 

One of the primary causes of the inefficient neutrino luminosity is the strong radial gradients of the 
entropy and leptonic fraction. Stable stratification of these quantities appears at the outer layer of PNSs 
and suppresses convective motions (Janka & Muller 1996). Convective regions are located deep inside of 
PNSs and surrounded by a convectively stable layer in which the surface fluxes of Vf, and ry are mainly built 
up (Buras et al. 2003). Thus the convection in PNS has little influence on the emission of neutrinos and is 
irrelevant for the supernova dynamics. 

On the other hand, MHD instabilities could be a seed of turbulence even in the stably stratified layer. 
Magnetic effects on the dynamics of supernovae have been examined actively since the importance of mag- 
netorotational instability (MRI) in supernova cores is pointed out (Akiyama et al. 2003). Two-dimensional 
MHD simulations of the collapse of supernova cores indicate that the shapes of shock waves and the neutrino 
sphere can be modified by the effects of magnetic fields (Kotake et al. 2004; Yamada & Sawai 2004; Takiwaki 
et al. 2004). Thompson et al. (2005) suggest that the turbulent viscosities sustained by the MRI convert 
the rotational energies of supernova core into the thermal energies and may affect the supernova dynamics. 
However, effects of magnetic phenomena inside of the PNSs have not been understood yet and still have 
many uncertainties. 

Therefore, in this paper we analyze the linear growth of MHD instabilities in the stably stratified layer 
of PNSs, and discuss the possibility of enhancement of neutrino luminosities. Especially, we focus on the 
instabilities caused by a strong differential rotation. We assume the toroidal component of the field is 
dominant in the supernova cores. Local stability of stellar objects with a rotation and toroidal magnetic 
fields have been investigated by many authors (Fricke 1969; Tayler 1973; Acheson 1978; Pitts & Tayler 
1985; Schmitt & Rosner 1986; Spruit 1999, 2002). Spruit (1999) concluded that the kink-type instability 
(Tayler instability) grows at first among various MHD instabilities in stellar radiative zone. However, most 
of the previous work assume a weak differential rotation because their motivation is for the interiors of main 
sequence stars like the Sun. The stability analysis for the cases of a strong differential rotation, which is 
expected in PNSs, has not been carried out sufficiently. In such cases, nonaxisymmetric magnetorotational 
instability (NMRI) could be dominant over the kink-type instability (Balbus & Hawley 1992). The nonlinear 
evolutions of the NMRI can initiate and sustain MHD turbulence and amplify magnetic fields, which may 
lead to the enhancement of neutrino luminosity via magnetoconvection, viscous heating, and the other MHD 
phenomena. 

The paper is organized as follows. In §2 we obtain the dispersion equation that determines the stability 
of PNSs. In §3 we discuss the general properties of the dispersion relation, and in §4 stability criteria and the 
maximum growth rate are derived by an analytical approach. In §5 we apply our results to the interiors of 
PNSs, and examine the linear growth of MHD instabilities. The conditions that the NMRI can be dominant 
over the kink-type instability and the nonlinear effects of MHD instabilities are discussed in §6. Finally, we 
summarize our main findings in §7. 
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2. THE DISPERSION RELATION 

We examine the stability of PNSs with a toroidal magnetic held using a linear perturbation theory. 
Acheson (1978) obtained a local dispersion equation for MHD instabilities in stellar objects with a toroidal 
held. We derive our dispersion equation following Acheson’s analysis. One of the original features of this 
work is that we consider the effects of leptonic gradients which can be important in PNSs. The size of 
differential rotation in PNSs could be much larger than that in other stellar objects (Villain et al. 2004). 
Therefore we focus on the cases of a strong differential rotation which are not well studied in previous work. 

The leptonization and cooling timescales of PNSs are much longer than the growth time of the instabili¬ 
ties considered in this paper. Thus the unperturbed state can be treated in a quasi-stationary approximation. 
The conductivity of plasma inside PNSs is so large that the decay timescale of magnetic fields is much longer 
than their cooling time. Then, the governing equations are the ideal MHD equations; 
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where p is the mass density, u is the fluid velocity, p is the pressure, B is the magnetic field, g is the gravity, 
and p, is the magnetic permeability. 

We use the cylindrical coordinates (r, (j>, z), and consider the Eulerian perturbations (denoted by a prefix 
S) with the WKB spatial and temporal dependence, S oc exp{i(krr + m</) + kzZ — at)}. Our analysis is purely 
local one, formally valid only in the immediate neighborhood of a particular location (r, z) and then only for 
perturbations of wavelength A in the r-z plane such that -k 2 ^. The azimuthal wavelength can be 

quite large, perhaps order of radius. 

Let us consider a PNS rotating with the angular velocity ft and toroidal magnetic helds B^. We assume 
these variables depend only on r and have a power-law dependence; f2(r) oc and B^{r) oc r~^. At the 
unperturbed state, the PNS is in the magneto-hydrostatic equilibrium, that is. 
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in the radial and vertical directions, respectively. 


When written out in component form and only the largest terms retained, the above set of equations 
becomes to linear order ; 
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Here we assume the local approximation {k ^ 1/r, m/r). We also adopt the Boussinesq approximation (e.g., 
Chandrasekhar 1961) because the Alfven speed and the rotation velocity of PNSs would be much smaller 
than the sound speed ~ 0.1c where c is the speed of light. In addition to these linearized equations, we use 
a relation by virtue of the local approximation, dp + = 0. This relation means that a fluctuation of 

the total pressure is negligible (Acheson 1978). 


To complete our set, we require another linearized equation. The chemical equilibrium is realized in 
PNSs, and thus the density is given by a function of pressure p, temperature T, and leptonic fraction Yi = 
(ue+n^)/n where Ue, n^, and n are the number densities of electrons, neutrinos, and baryons, respectively. In 
the Boussinesq approximation, the pressure perturbations are negligible because fluid elements are assumed 
to be moving slowly and in the pressure equilibrium with their surroundings. Then, the perturbations of 
density 5p and entropy per baryon 5S can be expressed in terms of the perturbations of temperature 5T and 
leptonic fraction dY): 


5p = -p{a^+S,6Yi'^ , 

Crp 

SS = msCp— + C6Yi , 


(13) 

(14) 


where tob is the baryon mass, a = — (9 In p/9 In T)pY,, ^ = — (91np/9y/)pT are the coefficients of thermal 
and chemical expansion, and Cp = {T/mB){dS/dT)pYi is the specific heat at constant pressure. The sign of 
(/ = {dS/dYi)pT determines whether the leptonic gradient is stabilizing or destabilizing. Here we neglect the 
effects of thermal and chemical diffusions for simplicity. The role of these diffusive processes on the growth of 
MHD instabilities is discussed later in § 6.1. The equations governing the thermal balance and the diffusion 
of leptonic fraction are written as 

f)Q 
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^+u-VYi=0, 

ot 


(16) 


(Miralles et al. 2004) and then the linearized equations are given by 
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Substituting equations (17) and (18) into equation (13), we obtain the linearized equation, 
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Eliminating perturbed quantities in eight linearized equations [eqs. (6) - (12) and (19)], we obtain the 
following fourth-order dispersion equation. 
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and 


are the square of a poloidal wavenumber, the square of the epicyclic frequency, the Alfen frequency, and 
the effective gravity vector. When we assume that isobaric and isochoric surfaces coincide, the dispersion 
equation (21) is simplified further by noting that 


g* X = 0 . 


(26) 


Then, defining = g*Wr'i’ with an analogous expression for N'^, the final form of the dispersion equation 
becomes, 
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where 


The quantities W and are components of the Brunt-Vaisala frequency N = {g* • V'l')^/^. As the leptonic 
gradient is considered here, the Brunt-Vaisala frequency is divided into the thermal buoyancy frequency Nt 
and the leptonic buoyancy frequency Nyi ; 


= N^ + , 
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where 


n4 
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If the leptonic gradient is ignored, the dispersion equation (27) is identical to that obtained by Acheson 
(1978) in the incompressible limit {cs —>■ oo). There are two physical branches to the dispersion equation: An 
internal gravity wave branch which is present in the absence of magnetic fields, and an Alfven wave branch 
which becomes unstable at sufficiently long azimuthal wavelengths when magnetic fields are present. 

Notice that the local analysis is irrelevant in the case that the growth rate of unstable modes is smaller 
than the differential rotation rate \dn/dlnr\ = qfl. In that case, this problem should be treated as an 
initial-value problem (Balbus & Hawley 1992; Kim & Ostriker 2000). This is because nonaxisymmetric 
disturbances cannot have a simple plane waveform in the presence of the shearing background on the wave 
crests (Goldreich & Lynden-Bell 1965). In this paper, we focus mainly on the NMRI for strong shear cases 
{q ~ 1) of which the growth rate is comparable to the angular velocity and use the local dispersion equation 
(27) for the stability analysis for simplicity. Kim & Ostriker (2000) examine the linear growth of the NMRI 
using both the local WKB analysis and direct integration of the linearized shearing-sheet equations taking 
account of the time-dependence of the radial wavenumber. They found that the instantenious growth rate in 
the direct integration is in good agreement with the WKB dispersion relation. The advantage of the simple 
WKB analysis is that we can derive analytical formula of the stability criteria and the growth rate of the 
instabilities. In fact, the stability conditions obtained by the WKB dispersion equation for a case of weak 
differential rotation (Acheson 1978) are quite useful and applied to the recent studies of the evolution of 
massive stars (Spruit 1999; Heger et al. 2005). 


3. GENERAL FEATURES OF THE DISPERSION RELATION 

In this section we solve the dispersion equation derived in the previous section and examine the char¬ 
acteristics of MHD instabilities in PNSs. There are four important parameters in the dispersion equation 
(27): the shear parameter q, the Brunt-Vaisala frequency N, the ratio between the angular velocity and the 
Alfven frequency id/ujA, and the power index of the field distribution s. The nature of MHD instabilities, 
such as the growth rate and characteristic wavelength, depends on these parameters, especially on the shear 
parameter q and the Brunt-Vaisala frequency N. 

The shear parameter is quite small in the Sun where the kink instability would be dominant among 
MHD instabilities (Spruit 1999). However the shear parameter can be of the order of unity in PNSs. In that 
case, the magnetorotational instability (MRI) can be important as is in rotationally supported accretion 
disks (e.g., Balbus & Hawley 1998). The growth rate of the MRI could be faster than that of the kink 
instability by many orders. MHD turbulence generated by the MRI can transport angular momentum inside 
of PNSs, which may affect their dynamical evolutions. 

In stably stratified regions, the Brunt-Vaisala frequency N is real and denotes the size of the negative 
buoyancy. When N is large, the growth of the axisymmetric MRI is suppressed significantly because the 
Brunt-Vaisala oscillation prevents the radial motion of fluids (Balbus & Hawley 1994). If the field is purely 
azimuthal, axisymmetiric modes are stable for the MRI and only nonaxisymmetric modes can be unstable. 
The nonaxisymmetric MRI (NMRI) in accretion disks has been examined (Balbus & Hawley 1992; Terquem 
& Papaloizou 1996; Kim & Ostriker 2000) and the operating mechanism is described in detail by Kim & 
Ostriker (2000). In this paper we investigate the stabilizing effects of N on the growth of the NMRI. 
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In what follows, we fix the values of the other two parameters. The ratio of angular velocity to the 
Alfven frequency and the field distribution are assumed to be VI/lua = 30 and s = — 2, respectively. These 
are about their typical values in rotating PNSs (see § 5). 


3.1. Effects of the Shear Parameter 

First we show the dependence of the growth rate of an unstable branch on the shear parameter q. The 
unstable growth rate is given by the imaginary part of w, i.e., 7 = 9(u;). In this subsection, the Brunt-Vaisala 
frequency is assumed to be = 0. Figure 1 is a contour plot of the growth rate 7 for a strong differential 
rotation case with q = 1.5. The growth rate is shown as a function of the azimuthal wavenumber m and 
the ratio of poloidal wavenumbers (3 = kr/kz- In this figure, the azimuthal wavenumbers are treated as real 
values instead of integer. As is expected, the NMRI modes can be seen for this case and the unstable growth 
rate is of the order of the angular velocity O. The maximum growth rate is 7 max ~ 0.7511 at /? = 0 and 
TO = 29. For a given to, the growth rate takes the maximum at /? = 0. The modes with a larger azimuthal 
wavenumber are stabilized with increasing (3. When kr ^ 0, the radial displacements of perturbations lead 
to undulations of the magnetic pressure, which acts as a stabilizing force for NMRI. Only when to is small 
there exists unstable modes with a larger (3^1, although the growth rate is much smaller than fl. 

Figure 2 shows the growth rate 7 as a function of the azimuthal wavenumber to for /? = 0, 0.5, 1, 2, 
and 4. Model parameters are the same as those in Figure 1. The shape of the dispersion relation is quite 
similar to that of the axisymmetric MRI (Balbus & Hawley 1991). For a given (3, unstable modes exist when 
the azimuthal wavenumber to is smaller than a critical value Wcrit and the growth rate takes the maximum 
at TOmax- These characteristic wavenumbers and the maximum growth rate decreases as l3 increases. When 
/3 = 0, the critical and the most unstable wavenumbers are TOcrit = 51 and TOmax = 29, respectively. 

It would be helpful to recall the axisymmetric MRI to understand the physical meanings of the char¬ 
acteristic values TOmax, Wcrit, and 7 max shown in Figures 1 and 2. The vertical wavenumber of the fastest 
growing mode fcmax and the critical wavenumber fccrit for the axisymmetric MRI are given by fcmax = 
a/— -I- Aq Q/2va and fccrit = -y/Sg ^/va, and the maximum growth rate is 7 max = gH/2 (Balbus & Hawley 
1991). When the vertical wavenumber k is replaced by m/r, we can estimate the azimuthal wavenumbers 
Wmax = a/— g^ -I- 4g Vl/2uja ~ 29 and TOcrit = -y/^g VL/uja ~ 51 for g = 1.5 and H/wa = 30. The maximum 
growth rate is given by 7 max = gH/2 = 0.75H for g = 1.5. These values are exactly consistent with those 
derived from our dispersion equation (27). 

Figure 3 shows the growth rate 7 as a function of the azimuthal wavenumber to for the cases with 
different shear parameters g = 0, 0.01, 0.1, 0.5, and 1.5. The Brunt-Vaisala frequency and the poloidal 
wavenumber are assumed to be V = 0 and (3 = 0. As seen from this figure, the growth rate of NMRI 
decreases as the shear parameter g decreases. When g ~ 1, the growth rate is comparable to the angular 
velocity. However, if g is less than unity, the maximum growth rate is almost proportional to g. The 
characteristic wavenumber TOmax becomes smaller as the shear is weaker. When g < 10“^, only to = 1 mode 
is unstable and the the growth rate is about 10 “^H. As is described above, the maximum growth rate of 
axisymmetric MRI is proportional to g, and thus it reaches to zero at g = 0. However nonaxisymmetric 
TO = 1 mode is still unstable even for the rigid rotation case (g = 0). This is quite different from the 
axisymmetric case and a unique feature of the NMRI. 



3.2. Kink-type (Tayler) instability 


The stability of the system with a toroidal magnetic field and a weak differential rotation is investigated 
by many authors (Tayler 1973; Acheson 1978; Schmitt & Rosner 1983; Pitts & Tayler 1986). They found 
that the kink-type instability (m =1) would be important in rigidly rotating stellar interior. The instability 
criterion for the m ^ Q mode depends on the radial distribution of the field, s = —dhiB^/dlnr (Spruit 
1999). The condition for the instability is given by 

(30) 


Obviously the kink mode (m = 1) 
pinch mode (m = 0) is given by 


is unstable when s = — 2. By comparison, the instability criterion for the 


s < -- 



(31) 


The pinch mode (m = 0) is usually stable in rigidly rotating stellar interiors if O loa- 


The growth rate of the kink-type instability is of the order of (Spruit 1999), which is about 10“^O 
for the case of VI/uja = 30. Thus the unstable mode in the limit of g ^ 0 shown in Figure 3 must be 
the kink-type instability. Notice that the NMRI and kink-type instability are on the same branch of the 
dispersion equation and the kink-type instability corresponds to an asymptotic solution of the NMRI in the 
limit of the rigid rotation (g ^ 0). The growth rate of the kink-type instability is much smaller than that of 
the NMRI in PNSs with a strong differential rotation, so that it would be anticipated that the NMRI would 
be dominant there. 


Hereafter, we describe the unstable modes with m = 1 as ’’kink modes” in distinction from ’’NMRI 
modes” of m > 2. The marginal shear rate is estimated by using typical growth rates of the NMRI modes 
(~ qfl) and the kink modes (~ w^/0). From the comparison of these growth rates, we can easily find that 
NMRI modes are dominant in the strong differential rotation cases q ^ (w^/fl)^, while the kink modes are 
dominant in the weak differential rotation cases q ^ {uja/^Y ■ 


3.3. Effects of the Brunt-Vaisala Frequency 

In this subsection, we investigate the effects of the Brunt-Vaisala frequency N on the NMRI and kink 
modes. The Brunt-Vaisala frequency N is real at stably stratified regions. The restoring force acts on 
displaced fluid elements along the direction of the entropy and leptonic gradients AI, which is usually toward 
the center of the star. Using a polar angle 9 = tan“^(r/z), the radial and vertical components of the 
Brunt-Vaisala frequency are written as 

Nr=Ns\rLe, N^ = Ncose. (32) 

Figure 4 shows the effects of N on the NMRI modes. The growth rate is shown by a gray contour on the 9-(3 
plane for the case of m = 29, which corresponds to the fastest growing mode of the NMRI. Since the size 
of the Brunt-Vaisala frequency in PNSs is highly uncertain (see § 5), we investigate a wide range of N/Q.. 
From top to bottom, Figure 4 shows the results for N/Q. = 0.1, 1, and 10. The shear parameter is assumed 
to be g = 1.5. 

The growth rate of the NMRI decreases with increasing N. The stabilizing effect by the Brunt-Vaisala 
oscillation depends on the angle between the poloidal wavevector fe and TV. When TV/0 = 0.1 (Fig. 4a), the 
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effect of N is weak, so that the growth rate is independent of 0. The mode with /3 = 0 is always the fastest 
and the maximum growth rate is 0.75fl. When N is comparable to (Fig. 4b), the growth rate of the modes 
with fc _L TV is slightly suppressed due to the restoring force of N. On the other hand, unstable modes with 
the poloidal wavevector parallel to N is unaffected by the Brunt-Vaisala oscillation and still have a growth 
rate of the order of the angular velocity O. If the Brunt-Vaisala frequency is much larger than 0 (Fig. 4c), 
the unstable modes can exist only when k || N. Near the equatorial plane {6 > 60°), the growth of the 
NMRI (to = 29) is completely suppressed. 

Next, we examine the effect of N on the kink modes. Figure 5 shows the dimensionless growth rate 
7/0 as a function of (3 and 9 for the case of to = 1. Notice that the growth rate of kink modes is much 
smaller than that of NMRI modes. When N is comparable or larger than fl, almost all modes except that 
k is parallel to N are suppressed, the same as in the case of NMRI modes. 

The difference between the NMRI and kink modes can be seen in the vicinity of the equatorial plane. 
The NMRI modes are suppressed around equatorial plane with increase of N/il (Fig. 4c). On the other 
hand, the kink modes can survive even near the equatorial plane (Fig. 5c). The direction of the entropy 
and leptonic gradients is almost horizontal when 6 ~ 90°, so that only a mode with a wavevector k || TV, or 
/3 = kr/kz 1, can grow against the Brunt-Vaisala oscillation. As seen in Figure 1, the kink mode can exist 
in the limit of /3 1, while the NMRI is suppressed completely. For the growth of the NMRI, the radial 

displacement of a perturbed fluid element is essential. However, the radial motion is stabilized near the 
equatorial plane and the NMRI cannot grow in such the situation. This would be the origin of the different 
behavior between two instabilities near the equatorial plane. 


4. ANALYTICAL TREATMENT 

In this section, we discuss the stability criteria and the maximum growth rate using an analytical 
approach. Since the dispersion equation (27) has a very complex form, it is difhcult to treat for what it 
is. Therefore we simplify it under realistic approximations. Our simplified dispersion equation captures 
the features of the original equation precisely, so that our analytical stability criteria would be useful for 
applications to the stellar interiors. 

The stability criteria for the cases of a weak differential rotation, q <C (w^/O)^, and a rigid rotation, 
<7 = 0, have already been investigated in terms of the interior of normal stars (Fricke 1969; Tayler 1973; 
Acheson 1978; Schmitt & Rosner 1983; Pitts & Tayler 1986; Spruit 1999, 2002). Here we consider a strong 
differential rotation case, q (w^/H)^, because it would be common in PNSs. In that case, NMRI modes 
are dominant and then we can use the high rotational frequency approximation assuming that the growth 
rate is comparable to the angular velocity and much larger than the Alfven frequency, w ~ H ^ uja- 


4.1. Stability Criteria 

The dispersion equation (27) can be rewrite as 

^ — [—2(1 — s)uj\ + + Ag] {up' — ttPuj'a^ — Auj\ {uo + toH)^ = 0 . (33) 

kz 

We consider a case of high rotational frequency, w ~ H ^ uja, and a large azimuthal wavenumber, to 1. 
Then the Alfven frequency in the second term can be negligible compared to the epicyclic and Brunt-Vaisala 
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frequencies. We can also approximate w + mO ~ in the third term. Using these approximations, we can 
obtain a simplified dispersion equation, 




+ fVg) = 0 . 


This is a quadratic equation of oj and the analytical solution is 

Ifc? " 


2 2 2 I 

OJ — 7Tl OJ ■\- 


2A:2 


1 1/2' 


+ N^±\ {k^ + Niy + 


(34) 


(35) 


The stability condition is given by > 0, that is, 

ft) ^2 +(^)ai+Ao>0, (36) 

where 

A 2 = + Nl , Ai = —2NrNz , Aq = — 2qQ? + m'^oj\ . (37) 

This inequality [eq. (36)] is satisfied if and only if (i) A 2 + Aq > 0 and (ii) B = A^ — 4 A 2 A 0 < 0 . By using 
these two conditions, we can derive two criteria for stability; 




(38) 


N'^oj\m'^ - 2qn^N^ > 0 . 


(39) 


These criteria are quadratic functions of the azimuthal wavenumber m. Thus, equations (38) and (39) are 
consistently satisfied for every azimuthal wavenumber in the case that zeroth order term of m is positive. 
Therefore the general stability criteria can be derived as 


^.2 n 

^ > 0 , 
omr 

(40) 


(41) 


These criteria have forms similar to the Solberg-Hpiland criteria but with gradients of angular velocity 
replacing the gradients of specific angular momentum. If the leptonic gradients are ignored, these criteria 
are exactly identical to those of the axisymmetric MRI (Balbus & Hawley 1995). 


4.2. Maximum Growth Rate 


The maximum growth rate of NMRI can be derived from equation (35). By differentiating partially 
with respect to m, the azimuthal wavenumber Wmax that provides an extremum is obtained as 

11/2 


1 


16 ^ 4 - + 


k 1 

Substituting rumax into equation (35), the maximum growth rate with respect to m is given by 


7n 


qft k. 




2 k 


1 Ni 


2q 


(42) 


(43) 
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We consider the characteristics of 7max in two limiting cases, = 0 and N ^ 0. When = 0, the 
unstable growth rate takes the maximum value, 


7 max = ) 


when kr/kz =0 and 


2uja 


\/-g2 + 4gr . 


(44) 


(45) 


These results are consistent with the characteristics of the original dispersion equation (27) discussed in 
§ 3.1. 

If A^ ^ 0, the ratio of the poloidal wavenumbers for the fastest growing mode /3max can be obtained by 
differentiating equation (43) with respect to (3. Then /3max should satisfy the following equation. 


- /3max^ f/^max + Pmax + “ “^qf^max f 1 “ ^ ■ 


(46) 


When N n, we can solve this equation approximately and the solution is given by /3max ~ Nr/Nz. This 
means that the poloidal wavevector of the fastest growing mode is parallel to the entropy and leptonic 
gradients {k || N). This is also consistent with the results of the original dispersion equation shown in 
Figure 4. 


5. APPLICATION TO PROTO-NEUTRON STARS 

In this section, we apply the simplified dispersion equation (34) to the interiors of PNSs and investigate 
the stability for the NMRI. In order to do that, we must know the amplitude of the physical parameters in 
the dispersion equation. Here we estimate physical quantities of PNSs based on the observations of pulsars 
and the numerical simulations of core-collapse supernovae. 

In the following we concentrate our discussions on the convectively stable regions in PNSs, which locates 
at the outer layer of PNSs (Janka & Muller 1996). If this region is magnetohydrodynamically unstable, 
neutrino luminosities would be amplified by various nonlinear magnetic processes (see § 6). The gain region 
behind the accretion shock is another site of MHD instabilities being important. Thompson et al. (2005) 
suggest that the viscous dissipation caused by MRI can enhance the heating in the gain region sufficiently 
enough to yield explosions. However, the gain region is convectively unstable. Convective motions may affect 
the growth of the MRI as is shown in convection-dominated accretion disks (Narayan et al. 2002). Although 
the linear and nonlinear interaction between the NMRI and convection would be a quite important issue, it 
is beyond the scope of this paper. 


5.1. Physical Quantities in Proto-Neutron Stars 

PNSs are expected to rotate differentially, as opposite to old neutron stars. Even when the iron core 
of a supernova progenitor is almost rigidly rotating (Heger et al. 2000; Heger et al. 2005), the collapse 
can generate a significant amount of differential rotation. Numerical studies of the rotating core-collapse 
indicate that the outer layers of PNSs would be rotating at the angular velocity 12 ~ 100 - 1000 sec“^ with 
a shear parameter g < 1 (Buras et al. 2003; Villain et al. 2004; Kotake et al. 2004). This angular velocity 
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is comparable to the observations of young isolated pulsars associated with supernova remnants (Marshall 
et al. 1998). Thus, we adopt a rotational profile of = 100 sec“^ and g ~ 1 as typical values in PNSs. 

MHD simulations of the core-collapse suggest that toroidal magnetic fields would be much stronger than 
the poloidal component (Akiyama et al. 2003; Kotake et al. 2004; Yamada & Sawai 2004; Takiwaki et al. 
2004; Ardeljan et al. 2005). This is because the toroidal component is generated by wrapping of poloidal 
fields by strong differential rotation during the core-collapse phase. On the other hand, observations of young 
isolated pulsars indicate that poloidal magnetic fields of neutron stars are typically of the order of G. 

Thus we expect that the toroidal fields at the interiors of PNSs could be 10^^ G. When we employ = 10^^ 
G, the Alfven frequency is 

UA ^ 3.0 f ^ - 3 ) sec“i . (47) 

G J VlO®-® cm/ g cm^^y 

The size of the Brunt-Vaisala oscillation in the stably stratified layer of PNSs is sensitive to microscopic 
physics such as the equation of state and neutrino opacity (Buras et al. 2003; Kotake et al. 2004; Thompson 
et al. 2005; Dessart et al. 2005). Therefore, we treat the Brunt-Vaisala frequency as a parameter. Because 
the Brunt-Vaisala frequency would be larger than the Alfven frequency in the stellar interiors, ^ 
the models with V/0 = 0.1 - 10 are examined here. We use our simplified dispersion equation (34) for the 
stability analysis. This equation does not include the parameter s, so that we can neglect the effects of field 
distribution in the following discussion. 


5.2. Stability of Proto-Neutron Stars 

The characteristics of the NMRI in stably stratified layers of PNSs are shown in Figure 6 . We assume 
the physical parameters q = 1, = 100 sec“^, and oja = 3.0 sec“^. For the Brunt-Vaisala frequency, we 

examine three cases with TV/fl = 0 . 1 , 1 , and 10 . 

Figure 6 a depicts the normalized azimuthal wavenumber of the fastest growing mode rrimax^A/^ ob¬ 
tained from equations (42) and (46). The characteristic wavenumber mmax depends on the polar angle 6 
and the Brunt-Vaisala frequency N. Figure 6 b shows the dimensionless maximum growth rate 7 max/fl as a 
function of the polar angle 9. The maximum growth rate 7 max is determined by /3max using equation (43). 
When the Brunt-Vaisala frequency is small N/n = 0.1, the azimuthal wavenumber rumax is independent of 
the polar angle 9. The maximum growth rate 7 max is constant at a value given by equation (44), which means 
that the NMRI is quite important in the entire regions. If the Brunt-Vaisala frequency N is comparable to 
the angular velocity N/D, = 1, the maximum growth rate decreases near the equatorial plane. However the 
decrease of the growth rate is at most by a factor of 2 at 0 ~ 90°. The azimuthal wavenumber is much larger 
than unity at any angle 9 so that the NMRI is still dominant over the kink-type instability. 

Even when N/Q is much larger than unity, all the regions are magnetohydrodynamically unstable. 
However the kink modes are dominant in the vicinity of the equator and thus the maximum growth rate 
decreases by many orders. Near the rotational axis, on the other hand, the azimuthal wavenumber Wmax 
and the maximum growth rate 7 max are unaffected by the Brunt-Vaisala frequency N. We define the critical 
angle dcrit at where the maximum growth rate becomes half of the value at 0 = 0. The critical angle is 
about 9 ~ 60° for the case of N/il = 10. We find that the critical angle is independent of the Brunt-Vaisala 
frequency if N/ft > 10. 

As a result, we can conclude that stably stratified layers of PNSs are always unstable to MHD instabili- 
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ties. For the cases of N/il < 1, the NMRI can grow at any angle 9. Even when N/fl is larger than unity, the 
polar regions is unstable to the NMRI. The growth rate is comparable to the angular velocity and much 
faster than that of the kink modes. The kink modes can be dominant only in the vicinity of the equator 
when the Brunt-Vaisala frequency exceeds the angular velocity. The typical growth time of NMRI is about 
100 msec. This is much shorter than the neutrino cooling time of PNSs (~ 10 sec), and thus the nonlinear 
growth of NMRI may affect the neutrino luminosities. 

We adopt Q./uia = 30 as a typical value in the discussions above. Recent numerical studies of rotating 
core-collapse suggests the possibility of more rapidly rotating PNSs up to ~ 1000 sec“^ (Fryer & Heger 
2000; Dimmelmeier et al. 2002; Thompson et al. 2005) and more strongly magnetized PNSs with ~ 10^^ 
G as an origin of magneters (Kotake et al. 2004; Yamada & Sawai 2004; Takiwaki et al. 2005; Sawai et al. 
2005). Therefore, it would be interesting to consider the rapidly rotating case, ^I/oja ^ 100, and strongly 
magnetized case, H/uja — 1 . 

The growth rate of the NMRI is proportional to the angular velocity. For rapidly rotating PNSs 
{il/ujA ^ 100), the NMRI can grow in shorter time scale than 10 msec. The characteristic length scale in the 
azimuthal direction, va/H, shifts to the shorter wavelength. In strongly magnetized PNSs {fl/uA — 1), on 
the other hand, the NMRI is stabilized by the strong magnetic tension. The critical azimuthal wavenumber 
for the growth of NMRI is given by ft/coA, and thus the NMRI modes would be suppressed completely 
if the field is very strong loa ~ 11. Only the kink modes can survive in such a situation because they are not 
affected by the magnetic tensions. 


6. DISCUSSION 
6.1. Effects of Diffusion Processes 


In the interior of PNSs, neutrino radiation plays an important role in the heat and lepton transports. 
The thermal and lepton diffusions caused by the neutrino radiation can reduce both the thermal and leptonic 
buoyant frequencies significantly. The coefficients of the thermal and chemical diffusivities in PNSs are given 
by 


10 ® 

M2O 


cm sec 


^chem — 10" 


/2 

^30 


II 

M20 


cm sec 


(48) 


where Yi/a is the electron fraction in units of 1/3, I 30 is a length scale in units of 30 km, T 4 is the temperature 
in units of 4 MeV, and /i 2 o is the electron chemical potential normalized by 20 MeV (Socrates et al. 2005). 


These diffusion processes must be considered if the wavelength of fluctuations is shorter than the critical 
length scale that can be estimated by equating the diffusion time and the period of the buoyant oscillations. 
The critical length scale of the thermal diffusion is given by At ~ ^/kt/Nt ~ 10® cm, where we assume 
Nt — 1000 sec“®. The critical length of the chemical diffusion may be a few times larger than that of the 
thermal one, because the leptonic buoyant frequency is smaller than Nt (Thompson & Murray 2000). Thus, 
our analysis is valid for the wavelength of the perturbations longer than the diffusion length At. 


However, the diffusive processes could affect the growth of the NMRI for the shorter wavelength A < At. 
The primary effect of the thermal and chemical diffusions is to reduce the restoring force caused by the stable 
stratification (Acheson 1978; Spruit 1999). Therefore, the maximum growth rate of NMRI increases and 
approaches to the case of iV/H ^ 0 when A <C At. Thus, the growth of these shorter wavelength modes may 
be dominant in the stably stratified layers of PNSs. 
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6.2. NMRI Modes vs. Kink Modes 

The conditions that the NMRI modes dominate over the kink modes can be derived from our simplified 
dispersion equation (34). The NMRI modes grow faster if the azimuthal wavenumber mmax is more than 
2. The shaded (white) area in Figure 7a represents the NMRI (kink) dominant region in the q-9 diagram. 
The boundary of these two areas is given by where mmax = 2 is fulfilled in equation (42). In this figure, we 
assume N/il = 10 and ^I/oja = 30 which are the typical values in PNSs. The NMRI modes are dominant at 
most of the regions if the shear parameter q is relatively large {q > 10“^). The kink modes become important 
only when the shear parameter is extremely smaller than unity. The boundary curve is little changed by the 
size of the Brunt-Vaisala frequency. 

For comparison, a similar picture for the solar radiative zone is shown by Figure 7b. In the solar 
radiative zone, the Brunt-Vaisala frequency is much larger than the angular velocity iV/fl ~ 10^ (Spruit 
2002). Here we assume the Alfven frequency VI/loa = 5 as the typical value in the solar radiative zone. The 
parameter region where the kink modes are dominant extends to the larger q regime because of the smaller 
ratio of H/wa- Only when the shear rate is about unity, the NMRI modes can be important for this case. 
According to the helioseismology, the solar radiative zone exhibits a nearly rigid rotation (<? <C 1). Thus, the 
kink modes are dominant and Tayler-Spruit dynamo would be amplifying the magnetic fields (Spruit 2002). 


6.3. Nonlinear Growth of the NMRI 

The NMRI can grow in the stably stratified layer of PNSs regardless of the size of the Brunt-Vaisala 
frequency. The nonlinear growth of the NMRI could initiate and sustain MHD turbulence and amplify 
magnetic fields. As is well known in the contexts of accretion disks, MHD turbulence leads to the angular 
momentum transport. Thus the rotational configurations within PNSs evolve toward a rigid rotation. If the 
PNSs will rotate rigidly, the neutrino sphere become oblate and the neutrino luminosies could be enhanced 
in the polar region (Kotake et al. 2003). Thus, the polar region would be heated by neutrinos preferentially. 

Turbulent viscosity generated by the NMRI can convert the rotational energies into the thermal energies, 
so that the neutrino luminosities from the PNSs may be increased (Thompson et al. 2005). Ramirez-Ruiz 
& Socrates (2005) suggest that the heating of materials through magnetic reconnection may also lead to 
the enhancement of neutrino luminosity. We would like to stress that the magnetoconvection and magnetic 
Rayleigh-Taylor instability of magnetic flux tube would assist the mixing of materials in the layers where is 
heretofore thought to be convectively stable (Parker 1974, 1975, 1979). 

Although these effects are absolutely a matter of speculation, it is sure that these effects are not con¬ 
sidered in the present scenarios of core-collapse supernovae. Therefore, we anticipate that these multiple 
magnetic effects can enhance the neutrino luminosity sufficiently to yield explosions in models that would 
otherwise fail. Because these magnetic phenomena are essentially three-dimensional, we must perform three- 
dimensional MHD simulations to reveal those effects quantitatively. This is our future work. 


7. SUMMARY 

In this paper, we perform the local linear analysis of the stably stratified layers in stellar objects taking 
account of toroidal magnetic fields and differential rotation. We derive the local dispersion equation and 
apply it to the interior of PNSs. Our findings are summarized as follows. 
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1. NMRI modes are dominant when the rate of differential rotation q is large and the growth rate is 

comparable to the angular velocity. In the limit of rigid rotation {q 0), only m = 1 mode is unstable, 
which corresponds to the kink instability. The growth rate of the kink mode is given by and typically 

much smaller than that of the NMRI. 

2. The growth of NMRI is suppressed by the Brunt-Vaisala oscillation. However, because the stabilizing 
force is along the direction of the entropy and leptonic gradients N, the unstable modes can exist if the 
wavevector is parallel to TV. Thus the kink modes is dominant over the NMRI modes in the vicinity of the 
equator when TV/H is large. 

3. Analytical formula of the stability criteria for NMRI and the maximum growth rate are derived 
from our simplified dispersion equation. The stability criteria are completely identical to those for the 
axisymmetric MRI if the leptonic gradient is ignored. The maximum growth is also agreed with that of the 
axisymmetric MRI for the case of TV = 0. 

4. NMRI can grow in the interiors of differentially rotating PNSs in spite of the size of the Brunt-Vaisala 
frequency. The suppression by the Brunt-Vaisala oscillation can be seen only near the equatorial plane where 
the kink mode is dominating. The typical growth time of NMRI is about 100 m sec, which is much shorter 
than the neutrino cooling time of PNSs. 
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Fig. 1.— The unstable growth rate 7 = 5(w) is plotted by the gray contour on the m-f3 plane, where m 
is the azimuthal wavenumber and (3 = kr/k^ denotes the ratio of the poloidal wavenumbers. The growth 
rate is normalized by the angular velocity f2. Three contour lines indicate y/fl = 0.7 (solid curve), 0.4 
(dashed curve), and 0.2 (dot-dashed curve). In this figure, we assume the shear parameter q = 1.5, the 
Brunt-Vaisala frequency iV = 0, the ratio of the angular velocity to the Alfven frequency ^IjuiA = 30, and 
the field distribution s = — 2 . 
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Fig. 2.— The unstable growth rate 7 /n as a function of the azimuthal wavenumber m. Each lines are 
labeled by the ratio of the poloidal wavenumbers (3. Model parameters are the same as those in Figure 1. 
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Fig. 3.— The dependence of the growth rate on the shear parameter q. The growth rate is shown as 
a function of the azimuthal wavenumber m for the cases oi q = 1.5, 0.5, 0.1, 0.01, and 0. The poloidal 
wavenumber is assumed to be /3 = 0. The other parameters are the same as in Figure 1. The growth rate 
decreases as the shear parameter decreases. Only the m = 1 mode can be unstable for the rigid rotation 
case {q = 0). 
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Fig. 4.— The contour plot of the growth rate of the NMRI modes (to = 29) for the case of (a) N/il = 0.1, 
(b) TV/n = 1.0, and (c) TV/fl = 10. The growth rate is shown on the 9-[3 diagram. The polar angle 9 is 
defined as 0 = tan“^(r/ 2 :) and measured with a unit of degree. The contour lines in each figure denote 
7/12 = 0.7 {dot-dashed curve), 0.4 {dashed curve), and 0.2 {dotted curve). Model parameters are the same as 
those in Figure 1 except for the Brunt-Vaisala frequency N. 
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Fig. 5.— The contour plot of the growth rate of the kink modes (m = 1) for the case of (a) N/Q = 0.1, (b) 
N/Q = 1.0, and (c) N/Q = 10. The growth rate is shown on the 6-13 diagram. The polar angle 9 is defined 
as 6* = tan“^(r/z) and measured with a unit of degree. The contour lines in each figure denote — 0.055 
{dot-dashed curve)^ 0.05 {dashed curve)^ and 0.04 {dotted curve). Model parameters are the same as those in 
Figure 1 except for the Brunt-Vaisala frequency N. 



















- 23 - 



6 (degree) 



6 (degree) 

Fig. 6.— (a) The azimuthal wavenumber and (b) the maximum growth rate of the fastest growing mode 
as a function of the polar angle 9 for the cases with N/il = 0.1, 1.0, and 10. These are obtained from 
our simplified dispersion equation (34). We choose the model parameters as q = 1.0, ft = 100 sec“^, and 
UJA = ^ sec“^. 
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Fig. 7.— (a) The NMRI and kink dominant regimes in the diagram of the shear parameter q and the polar 
angle 9 for the case of the PNSs. The border of the two regions is where the maximum azimuthal wavenumber 
is rumax = 2 in equation (42). Model parameters for this case are N = lOfl and 11 = SOwa- (b) The same 
plot for the case of the solar radiative zone. Model parameters are N = lO^fl and fl = ^uja- 





